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Abstract. Whereas the short time behaviour of an unstable quantum mechanical system is well 

understood from its theoretical as well as experimental side, the long time tail of the very same 

systems has neither been measured experimentally nor is there a theoretical agreement on how 

to handle it. We suggest a possible way out of this unsatisfactory state of art. Theoretically we 

suggest that the correct spectral function entering the Fock-Krylov method to calculate the survival 

amplitude is proportional to the density of states of a resonance. The latter is essentially the energy 

derivative of a phase shift. As a bonus, we can connect the survival probability to scattering data 

via the phase shift. The method then not only establishes the spectral function, but is per se a semi- 

T^ ' empirical method to extract the large time behaviour from scattering data. 
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A SHORT HISTORICAL TALE OF QUANTUM DECAY 

T— I ■ 
> 
(^ I This year, in 2008, we celebrate the fiftieth anniversary of a famous quantum mechan- 

O . ical result which has spawned numerous papers on the subject. In 1958 Leonid Khalfin 

proved that in the quantum world, the exponential decay law P{t) = e^^', is only an 
approximation [1 J . The quantum mechanical survival probability, without any approxi- 

(^ • mation is, 

g: P{t) = \A{t)\^ = me-'"'\yf)\\ (1) 

It can be calculated at short times to give, P{t) ~ 1 — (Aip//)?^, which in turn is connected 

k> ! to -^1^=0 = ^^ ^^^j ' \[=o 7^ 0. A direct deviation from the exponential decay law 

^ I at short times has been experimentally verified [2]. It is therefore correct to say that 

we have a very good understanding of the first few moments in the life of an unstable 
quantum state. This is not the case for the large time tail as not only has it never been 
experimentally seen (in spite of efforts to detect it [3]), but its theoretical foundation 
seems to 'enjoy' different treatments i^. Our idea presented in this talk is to use a 
semi-empirical method (experimental data in a theoretical formula) instead of a direct 
evidence for the large time behavior of the survival probability. We will show that this 
is also closely related to pinning down a more exact and reliable theoretical framework. 
To achieve the goals we need two different time concepts in Quantum Mechanics: (i) 
Time as a parameter and (ii) Time as an observable. The first is clearly the variable time 
t which appears in the Schrodinger equation and the survival probability; the second has 



to do with quantum mechanical observables[l| such as delay-time, dwell-time or sojourn- 
time, Larmor-time, traversal-time etc. which were constructed to answer questions about 
the time spent in a region or the quantum collision time. In calculating the parametric 
time dependence, the 'observable time' will be necessary. 

LARGE TIME BEHAVIOUR OF THE SURVIVAL AMPLITUDE 

The Fock-Krylov Method 

The Fock-Krylov method [5] is a suitable theoretical framework to study the large 
time behaviour of unstable systems. It relies on basic quantum mechanical results and 
therefore is, up to a point which we will discuss later, model-independent. We first 
observe that an unstable state |^) cannot be an eigenstate to the energy i.e. //|^) 7^ 
£■1^). Otherwise the survival amplitude A (?) and the survival probability P{t) would 
come out trivially to be, A{t) = (^|e^'^'|^) = e^^^' and P{t) = 1. Hence, assuming a 
continuum, H\E) = E\E), {E'\E) = d{E' — E),weaie entitled to expand 



1^)- / dEa{E)\E) (2) 

JSpect(H) 

where, 

p{E)^^-^^-mf)\'=\am\ (3) 

is a probability density (and as such positive-definite) to find the states with energy E 
in the resonance. This distribution is also known as the spectral function. We can now 
calculate the survival amplitude to obtain 

A{t)= f dEp{E)e^'^' = r p{E)e-'^' (4) 

JSpect(H) ''Eth 

which turns out to be a Fourier transform of the spectral function. E^h is the sum of the 
masses of the decay products. The success of this general method hinges on the right 
choice of p. There is no general agreement in the literature on what this function should 
be (in the next section we make a claim about the correct choice of p), but a general 
parameterization looks like, p{E) = (Threshold) x (Pole) x (Form — factor) , i.e., 

p{E) = {E-ErhyxP{E)xFiE). (5) 

Some comments about p are in order: (i) P{E) has a simple pole at zr = Er — ITr/I 
which leads to the exponential decay law. More poles in the fourth quadrant of the 
complex E-plane would modify even the exponential part of the decay, (ii) F{E) has 
no threshold and no pole. It is a smooth function which should go to zero for large 



^ Due to Pauli theorem there is no time operator as a conjugate variable to energy, but time-delay operators 
can exist nonetheless. 



E. (iii) Large times t correspond in the Fourier transform to small E (in agreement 
with the time-energy uncertainty relation). Hence, the large time behaviour is due to the 
choice of 7 which is often controversial, (iv) The transition region between exponential 
and non-exponential is partly due to the choice of the form- factor F . (v) The choice 
of the spectral function is not unique in the literature: e.g. often 7=0 and f{E) = 1. 
The problem regarding the right choice of the spectral function will be discussed in the 
next section. But even without the explicit knowledge of p, one can extract valuable 
information from the Fock-Krylov method. We choose the closed path in the complex 

£'-plane: Cr = C3 +C^+C^ , starting from zero along the real axis (Cg^) attaching to 

1/4 

it a quarter of a circle with radius R (C^ ) in the clockwise direction and completing the 
path by going upward the imaginary axis up to zero (C3). Hence, 

where the dots indicate the integrand from equation (|4]) with an argument shifted by 
Eih since we start from 0. The integral we wish to calculate is along the real axis. / is 
calculated by using the residue theorem with the pole at zr leading to the exponential 
decay law. We assume that forR -^ 00 the integral along the arc goes to zero. Thus, 

A{t) = AE{t)+AtT{t) 
AE{t) = 2niPizR)F{zR){zR-E,hre'''''e-^''/^ = aE{t)e-^''/^ (7) 

with P{z) = lim^^j.^ P{z) {z - zr) and 

POO 

^LT{t) — (phase) X / dxP{~ix + Eth)F {—ix -]- Eth)x^ e^^' 

Jo 

~ (phsisc) xr{Y+l)P{ETh)F{Erh)xj^ = aLT-^ (8) 

The result agrees with [g]. This is how nature slows down the exponential decay. 



The transition region and critical times 

One can approximately estimate the transition time from the exponential to the power 
law behaviour by setting |a£;|e~^'"^'o/^ ~ I^Lrl^qr or alternatively by determining the 

zeros of the function 

a)(^o)=ln^f^V +(7+l)ln^o-^o, ^0 = ^ (9) 

Note that, strictly speaking, there are three regions (see Fig. 1) and two critical times: the 
time at the transition from the exponential to the oscillatory region and from the oscilla- 
tory to the power law. The condition ^ can also have two zeros out of which normally 



- Exponential 




20 40 60 80 100 

rt 



FIGURE 1. Characteristic behaviour of a survival probability. 



the second one is the right indicator of critical time. Consider a narrow resonance model 

{E-ErY+(Tr/2)^ ■ 



with the spectral function a simple Breit-Wigner, p{E) o^ {Er — E, 
Then defining £ = Er —£■//, and choosing 7=1/2 



th) 



|«£| 



e-'/^^[l + e^/{rR/2)^] 
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|aLr|(r^/2)3/2-^ ,.,./,.„/^, , . (10) 

It is evident that smaller the ratio e/Tr, smaller is the critical timeH. Hence for narrow 
resonances the best candidates to find the non-exponential long time behaviour are 
threshold resonances. 



THE CONNECTION TO SCATTERING DATA 



Beth and Uhlenbeck [7], while calculating virial coefficients 5,C in the equation of an 
ideal gas: pV = RT 1 + f^ + ^ H , found that the difference between the density of 

states (of scattered particles) with interaction dni{E)/dE and without drii ' {E)/dE is, 

dn _ dnijE) dnf\E) _ 2/+ 1 d5i{E) 

dE dE dE K dE 

In a resonant scattering [8] this is the density of states of a resonance (in terms of the 
decay products). For instance T = {Tr/2)/{Er — E — iTr/I), gives, 

^ = Ylll (12) 

dE {Er-eY + tI/a- 

This offers a semi-empirical method to examine the large time behaviour of of unstable 
systems (resonances) directly from data since the spectral function is. 



P{E) 



JProb 
dE 



1 dn 
ndE 



(13) 



This is an approximate estimate since there are two critical times characterizing the transition region 



Note that this also uniquely fixes the spectral function, at least in the vicinity of a res- 
onance. This interpretation works well for all I- values except for the 5-wave, because 
in this case ^ »« ^ ,,^2 and we encounter a threshold singularity. This unreason- 
able singular behaviour of the density of states can be remedied without changing the 
interpretation (see next section). Otherwise for / > we have ^°^ [E —EthY^^'^ and 
therefore 7= / — 1/2, i.e., 7 is fixed. An explicit example of /'(f) for narrow resonances 
a + a ^ ^Be{l = 2) — > a + a has been calculated in [|9i] from scattering data given in 
form of the phase shift. 



EXPLICIT EXAMPLES 

The Fock-Krylov method is not limited to narrow resonances. We can put forth the 
question if there exist new features in the survival probability for broad resonances? 
For narrow resonances we can use a Breit-Wigner model for the amplitude to calculate 
the critical time and obtain some relevant results. No general parameterization of the 
transition amplitude for broad resonances exists. We will therefore choose the most 
prominent example of the o for which some parameterizations are available. The new 
features which can appear here are e.g., sub-threshold zeros (called Adler's zeros) in the 
amplitude and hence also in the form-factor in the density of states. Besides this, the o 
is an s-wave resonance. As we already mentioned the s-wave density of states has to be 
modified from the time delay to the dwell time delay. 

— = dwell - time = 2— y^- (14) 



which is the relativistic version of the expression found in [llOll . This singularity-free 
expression is also a density of states as shown in [1 IJ. There exist different parameteri- 
zations of the amplitude for TtTt ^ o ^ TtTZ. Here we have opted for the following one 



from [|12|1 



^^W (15) 



M^-s-iMY{s)' 

where the energy dependent width can be found in [|l2n . It has the structure: threshold x 
Adler's zeros x form-factor. Figure 1 displays the critical times (ft)((^o) from equation 
(|9l)) for very narrow resonances in nuclear physics (^Be) in contrast to the same exercise 
done for the broad o. The critical lifetimes of these narrow resonances are indeed large 
(~ 30 and 70 in terms of lifetimes). However, in the case of a broad resonance, i.e a, 
we can conclude the following, (i) The short transition time which is usually considered 
an artifact of the approximation (recall that in reality there are two critical times) is the 
biggest we find (ca. one lifetime). The large transition time which usually is the correct 
time scale for transition is the smallest we find (eight lifetimes). Therefore, we think that 
the oscillatory transition region could here be of importance and the real transition time 
could be somewhere between the two we find now. This behaviour would be similar 



in other broad resonance systems such as the rj-mesic nuclei [|13l] . The relatively small 
critical time and the importance of the oscillatory region makes the study of the time 
evolution of the o an interesting undertaking which we plan to continue in future. 
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FIGURE 2. Critical times for narrow resonances and O as an example of a broad resonance. 
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